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ABSTRACT 

In a low-metallicity gas, rapid cooling by dust thermal emission is considered to induce cloud 
fragmentation and play a vital role in the formation of low- mass stars (< 1 M ) in metal-poor 
environments. We investigate how the growth of dust grains through accretion of heavy elements in the 
gas phase onto grain surfaces alters the thermal evolution and fragmentation properties of a collapsing 
gas cloud. We calculate directly grain growth and dust emission cooling in a self-consistent manner. 
We show that MgSi03 grains grow sufficiently at gas densities hh = 10 10 , 10 12 , and 10 14 cm~ 3 for 
metallicities Z = 10~ 4 , 1CU 5 , and 10~ 6 Zq, respectively, where the cooling of the collapsing gas cloud 
is enhanced. The condition for efficient dust cooling is insensitive to the initial condensation factor 
of pre-existing grains within the realistic range of 0.001-0.1, but sensitive to metallicity. The critical 
metallicity is Z cv a ~ 10~ 5 ' 5 Zq for the initial grain radius 7"MgSiO 3 ,0 iS 0.01 /mi and Z CI i t ~ 10~ 4 ' 5 Zq 
for rM g sio 3 ,o ^ 0.1 /im. The formation of a recently discovered low- mass star with extremely low 
metallicity (< 4.5 x 10 -5 Zq) could have been triggered by grain growth. 

Subject headings: dust, extinction — galaxies: evolution — ISM: abundances — stars: formation — 
stars: low-mass — stars: Population II 



1. INTRODUCTION 



Dust grains in the early universe are considered to be 
crucial for low-mass star formation in the low-metallicity 
gas. 3 The efficient cooling by dust thermal emission 
can trigger the fragmentation of the gas at densities 
n H ~ 10 13 -10 15 cm~ 3 and gas temperature T g — 1000 K, 
where th e corresponding Jeans mass is less than on e so- 
lar mass ([Schneider et al.l 120031 : lOmukai et al.ll2005l ). In 
the early universe, dust grains are considered to be pre- 
dominantly supplied by core-collapse s upernovae; grains 
are formed in the supern ova eject a (|Todini fe Ferraral 
120011: iNozawa et ail 12003(1 but a part of them are de- 
stroyed by the rever s e shock penetrating into th e ejecta 
(jNozawa et alJl2007t iBianchi fe Schneider! l2007f ). Thus, 
the depletion factor of dust, /d ep (the mass ratio of dust 
to metals), is determined by the balance between the 
fo rmation and des t ructio n of grains. 

ISchneider et al.1 (J2012) followed the thermal evolution 
of gas clouds adopting the depletion factors for their 
models of dust formation and destruction in super- 
novae. For metallicity Z, the absolute amount of dust 
is related to th e dust -to-gas mass ratio, V = fd ep Z. 
ISchneider et al.l (|201 2f ) showed that the condition of 
cloud fragmentation owing to rapid dust cooling is de- 
scribed as the minimal dust-to-gas mass ratio, £> cr it, un- 
der the assumption that the dust-to-gas mass ratio is 
consta nt in the collapsing cloud. However. INozawa et al.l 
(2012) recently argued that the dust abundance can be 
increased by the growth of dust grains through accretion 
of heavy elements in the gas phase onto grain surfaces. 
They found that the condensation efficiencies (the num- 
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ber fraction of an element depleted onto grains) of Si- and 
Fe-grains increase up to unity at gas densities tt-h = 10 10 - 
10 14 cm -3 for elemental abundances [Si, Fe/H] = —5 to 
—3. However, they did not follow the thermal evolution 
of the gas clouds to determine whether the conditions for 
fragmentation are met in the clouds. In this Letter, we 
investigate the effects of grain growth on the thermal evo- 
lution of star-forming clouds by calculating grain growth, 
formation of hydrogen molecules on the grain surfaces, 
and dust cooling self-consistently. Then, by investigating 
possible range of metallicities and initial condensation 
efficiency factors, we can determine the critical metallic- 
ity, the minimal metallicity required for the formation of 
low-mass fragments. 




Fig. 1. — Planck mean opacities per unit dust mass for MgSi03 
(solid), Mg2Si04 (dashed), Fe304 (dotted), Fe (dot-dashed), and 
Si (dot-dot-dashed) as a function of dust temperature for a grain 
radius 0.01 fim. The values are insensitive to the grain radius 
except for Fe and Fe304. 
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Fig. 2. — Condensation efficiency (top) and the radius of MgSi03 
grains (middle), and gas temperature (bottom) as a function of 
hydrogen number density in gas clouds with metallicities 10 -4 Zq 
(solid), 1CT 5 Z Q (dotted), and 10 -6 Z Q (dot-dashed) for the initial 
condensation factor 0.1 and the initial grain radius 0.01 (im. With 
(red) and without (blue) grain growth. Filled circles mark the 
states where the fragmentation conditions are met (see the text). 

2. NUMERICAL METHOD 

We explore the thermal evolution of a collapsing gas 
cloud with metallicities of Z = lCT 6 -lCr 3 Z Q . The evo- 
lutions of density p and gas temperature T g of the cloud 
core are followed by a one-zone model. Our calculations 
include chemical reactions of gas species and radiative 
transfer for line emissions from the g as, and continuum 
emissions from the gas and dust (see lChiaki et al.ll2013l 
for details). 

In our calculations where the growth of dust grains 
is included, the mass density of dust pd increases with 
time. Thus, we calculate the continuum optical depth as 
follows; 

Tcont = (KgPg + KdPd) iBE, (1) 

where n g and denote the Planck mean opacities of gas 
and dust, respectively, 4 p g denotes the mass densities of 
gas (p = p g + pd) and ^be is the Bonnor-Ebert length 
as the typical size of the cloud. Planck-mean opacities 
of dust Kd are calcula ted using optical con stants from 
references tabulated in iNozawa et al.l (|2008T ) and shown 
in Figure [1] for grain species considered in this study. We 
then obtain the escap e fraction of co ntinuum emission as 
cont = mm{l, r-jt} (|Omukail l20M). 

In this study, we primarily consider the growth 
of MgSiC>3 grains; Mg-silicate is one of the ma- 
j or grain species eje cted b y core-collapse supe rnovae 
(|Bianchi fc Schneiderl 120071 INozawa et al.l 12007ft . We 
also investigate the growth of Mg 2 Si04 grains and its 

4 We should note that here the definitio n of dus t Pla nck mean 
opacity k<j is different from that of K gr in IQmukail 120001') in that 
the former is defined as the absorption cross section per unit dust 
mass, while the latter is defined as the absorption cross section per 
unit gas and dust mass. 



effects on the thermal evolution of the gas clouds (see 
Section Q] below) . We assume that the growth of these 
multi-element grains is regulated by the least abundant 
gaseous species amon g reactants, which we ca ll the key 
element hereafter (see iHirashita fc Kuo Il20ilft . 

Using the concept of the key element, we calcu- 
late grain gr o wth by following the prescription in 
INozawa et al.l (|2012f ). On the assumption that pre- 
existing dust grains are spherical with a single initial 
radius o, the evolution of grain radius, r^, is described 

by 

§-(t-S.)(^)">M. M 

where Si denotes the sticking probability of the gaseous 
atoms incident onto grain surfaces, a^o are the average 
radius of a monomer molecule in the dust phase per key 
element, and rrij and n^ as are the mass and the num- 
ber density of key element j, respectively. In equation 
@, the evaporation term can be neglected. When the 
cloud fragmentation is triggered through efficient cool- 
ing by dust, the temperature of dust increases only up 
to several hundred kelvin, well below the sublimation 
temperature of MgSiOs (^1000 K under the condition 
considered here). Therefore, dust evaporation does not 
affect the fragmentation condition in the situations we 
investigate. 

The time evolution of condensation efficiency /y, de- 
fined as the number fraction of the key element j 
locked in grain species i, is described by fij(t) = 
/ij,o[^i(i)/^i,o] 3 j where /y^o represents the initial conden- 
sation factor. Then, the mass density of the grains is 
given by Pd,i(t) = pX^Ajfijpi, where Xh is the mass 
fraction of hydrogen nuclei (set 0.76 here), and pi is the 
molecular weight of grains per key element. 

We set the initial density and temperature of a gas 
cloud to be n H ,o = 3000~ cm~ 3 and T g:0 = 200 K, 
respectively. 5 The initial number fractions of chemi- 
cal species relative to hydrogen nuclei are A c = 10~ 5 , 
A B = 2.87 x 10~ 5 , An 2 = IQ- 3 , A KU = 10~ 6 , irrespec- 
tive of the metallicity Z. The number abundance of an 
element j is given as Aj = A^q(Z j 'Zq), where the solar 
numb er abundanc e s of he avy elements, A, ©, are taken 
from lCaffau'eTaLl (|2011a[ ). 

The average radius of a n atom per key elem ent in the 
dust phase is taken from (|Nozawa et al.ll2~003f ). and the 
sticking coefficient is assumed to be unity. For the scaled 
solar elemental abundances (^Mg/^Si > l)j the key ele- 
ment of MgSi03 is Si, whose initial condensation factor 
is taken as /M g SiO 3 ,Si,0 = 0.001, 0.01, and 0.1. As for 
the initial radius of pre-existing grains, we consider four 
cases of rM g SiO 3 ,0 = 0.001, 0.01, 0.1, and 1 pm. Dust 
temperature is calculated from the balance equation be- 
tween heating due to collisions with the gas and cooling 
by dust thermal emission. Also, the formation rate of hy- 
drogen molecules on grain surfaces is calculated by con- 
si dering impact of hydrog en atoms onto grain surfaces as 
in lSchneider et alJ (|2006l ). 

5 The initial condition of the gas cloud is taken from the values 
when a supernova shell becomes gravitationally un stable in our 
one-d imensional calculations of a supernova remnant llChiaki et alJ 
120131 1. We have confirmed that the evolutions of collapsing gas 
clouds at rifj > 10 s cm -3 are insensitive to these initial values. 
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Fig. 3. — Fragmentation properties of the collapsing gas clouds for each model on Z-/MgSiO 3 ,Si,0 (left) and Z-?"MgSiO 3 ,0 (right) planes. 
For models marked by double circles, the fragmentation conditions are met both with and without grain growth. While, for models marked 
by filled circles, the conditions are met only with grain growth. For models marked with open circles, the conditions are not met both with 
and without gra i n gro wth. In left panel, the dashed line shows the metallicity corresponding to the critical dust-to-gas mass ratio given by 
Schneider ct al. (2012) (sec the text) and the solid line shows the critical metallicities obtained by detailed calculations. 



3. RESULTS 

Figure [2] shows the evolutions of the condensation 
efficiency /MgSi0 3 ,Si(i)> the radius r Mg siO 3 (0; and the 
gas temperature T g for the initial condensation factor 
/M g SiO 3 ,Si,0 = 0.1 and the initial grain radius r Mg siO 3 ,0 
( ! . ( ) 1 /im in gas clouds with metallicities 10 _6 -10 -4 Z Q . 
For the models with grain growth, the condensation ef- 
ficiency (mass density of grains) increases at nn = 10 10 , 
10 12 , and 10 14 cm" 3 ' for Z = 1CT 4 , 10" 5 , and 1(T 6 Z Q , 
respectively. It eventually converges to unity regardless 
of metallicities. This is because the timescale charac- 
terizing grain growth eventually becomes smaller than 
the timescale chara cterizing the free-fall collapse (see 
iNozawa et aTll2012f ). The grain radius also reaches to 

a single value r Mg si0 3 ,o(l//M g Si0 3 ,Si,o) 1/3 - Since the in- 
crease in grain mass density causes enhanced dust cool- 
ing, the models with grain growth lead to lower gas tem- 
peratures than the models in which the condensation ef- 
ficiency is constant (i.e., without grain growth). 

The circles in Figure [2] depict the phase (on each tra- 
jectory) where the cloud is expected to fragment (see 
ISchneider fc Qmukail I2010L for the detailed description 
of the criterion for fragmentation). Interestingly, for 
Z = 10~ 5 Zq, the trajectory reaches this phase with 
grain growth, whereas, without grain growth, the trajec- 
tory does not. We therefore conclude that grain growth 
can reduce the metallicity (Si abundance) above which 
low-mass fragments can form. 

3.1. Critical Metal Abundance 

Figure [3] shows fragmentation properties for each 
model we investigate. The double and filled circles de- 
pict the models for which the fragmentation conditions 
are met when we consider grain growth, and double cir- 
cles depict the models for which the conditions are met 
without grain growth. From these results, we can dis- 
cuss the critical condition for the formation of low-mass 
fragments. 

First, we make the comparison of our results with- 
out grain growth with pre vious works. As the m inimal 
condition of fragmentation, Schneide r et al.l ()2012l ) found 



the critical dust-to-gas mass ratio X> cr it = 4.4 x 10 -9 ( 
= fdepZ), above which the dust-induced fragmentation 
occurs. In the dust model they used, the grains have 
size distribution functions with the peak at 0.01 /im. 
Thus, we can compare their result with our models for 
r M g SiO 3 .0 = 0-01 A*m. To apply this, we should convert 
/M g SiO 3 ,Si,0 into /dep as 

, Pd,M g Si0 3 ^H^Si^M g Si0 3 t /Q \ 

Jdcp = — 7} /M g SiO 3 ,Si,0, V>) 

Pmctal & 

where p mo tai is the mass density of metal in both gas 
and dust phase and /iM g Si0 3 = 100. With th e so- 
lar metallicity Z Q = 0.0153 (jCaffau et al.l 1201 lbl ). we 
obtain / dcp = 0.16/ Mg siO 3 ,Si,0- In Figure El we plot 
Z = 2>crit/(0.16/ Mg si0 3 ,Si,o) by the dashed line. Our 
results are co nsistent with the critica l dust-to-gas mass 
ratio found by ISchneider et ail (|2012l ). 

Then, we discuss how the critical condition is changed 
if we consider grain growth. The solid line in Figure 
[3] shows the critical line above which the fragmentation 
condition is met with grain growth. Comparison of the 
two critical lines (solid and dashed line) shows that lower 
metallicity is generally allowed for low-mass fragment to 
form when we consider grain growth. Also, the solid 
line in Figure [3] shows that the critical metallicity is 
insensitive to the initial condensation factor. This is 
because even if the initial condensation factor is small 
(/M g SiO 3 ,Si,0 ~ 0.001), the grains grow up to unity be- 
fore the dust emission cooling becomes effective unless 
the metallicity is too small. 

Instead of the condensation efficiency (dust-to-gas 
mass ratio), the critical condition is determined by the 
metallicity and dust size. The estimated critical metal- 
licity is -2/crit — 10~ 5 ' 5 Z Q for the initial radius of MgSi0 3 
grains r Mg SiO 3 ,0 < 0.01 /im, and Z crit ~ 10~ 4 5 Z Q for 
/M g SiO 3 ,Si,0 > 0.1 /im. For larger initial radius, the gas 
density at which grains grow is higher because the to- 
tal surface area of dust particles, which determines the 
growth rate, is proportional to (for a fixed mass den- 
sity of dust). The right panel of Figure [3] shows that for 
r M g SiO 3 ,0 = 0.001 and Z = 10~ 3 Z Q , the fragmentation 
conditions are not met without grain growth, since heat- 
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Fig. 4. — Fragmentation conditions of collapsing gas clouds on 
/MgSiO 3 ,Mg,0 _r MgSiO 3 ,0 plane for the abundance of heavy elements 
of SDSS J102915+172927. Symbols represents the same meaning 
as Figure [3] 

ing due to H2 formation compensates for cooling by dust 
emission. 

3.2. Application to SDSS J '102915+172927 

We perform the calculations of grain growth and evolu- 
tion of collapsing gas, adopting the abundance of heavy 
elements for the star SDSS J102915+172927: A c = 
1.58 x 1CT 8 , A = 2.51 x 1CT 8 , A Mg = 8.91 x 1CT 10 , 
A Si = 1.78 x 10~ 9 , and A Fc = 3.39 x 1CT 10 , correspond- 
ing to the metallicity 4.5 x 10~ 5 Z Q . Note that the 
abundance of carbon is an upper limit because no strong 
carbon enhancement is evident in the stellar spectrum. 
Also, they estimate the oxygen abundan ce, assuming 
the typ ical excess of the a-element oxygen (jCaffau et al.1 
l2011al) . Given the abundance of the gas, we aim at con- 
straining the initial condensation factor and the initial 
grain radius necessary for triggering the formation of the 
low-mass star. 

Figure [4] shows the results of calculations with the 
growth of MgSi03 grains. In the present case, the 
key element is Mg because of ^Mg/^Si < 1- In- 
terestingly, for most of the cases considered here (all 
cases except for /MgSiO 3 ,Mg,0 = 0.01 and 0.001 with 
nvigSiO 3 ,0 = 1 Mm), the fragmentation conditions are 
met. Thus, we argue that the parent cloud of SDSS 
J102915+172927 could have been enriched with rela- 
tively small grains (rM g SiO 3 .0 ^ 0.1 /im) or other- 
wise with relatively high condensation factor of dust 

(/MgSiO 3 ,Mg,0 ^ 0.1). 

It is worth mentioning that the above results are 
approximately consistent with those obtained with the 
scaled solar composition (Figure although the abun- 
dance pattern is different. The results for rM g siO 3 ,0 = 
0.01 /im in Figure H] lie between those with metallicity 
log(Z/Z Q ) = —5 and —4 in the left panel of Figure El 
and the results for /MgSiO 3 ,Mg,0 = 0.1 lie between those 
with metallicity log(Z/Z@) = —5 and —4 in the right 
panel of Figure El 

4. DISCUSSION 

In our calculations for the scaled solar abundances, the 
least abundant element, silicon, regulates the growth of 
MgSiOa grains. Though we have assumed that Si atoms 



are directly accreted onto the grains, a part of or all 
the Si atoms might be in the form of SiO molecules in 
dense clouds. We also examine a model in which all Si 
atoms are in the form of SiO molecules initially. As a 
result, the density at which grains grow rapidly and the 
resulting critical metallicity are not changed from the 
results of models in which we consider Si atoms as the 
key element. 

We also examine the growth of other grain species and 
the effect of each species on the thermal evolution of 
clouds. We find that, for the scaled solar composition 
(Asi/Af c = 1), Fe304 grains grow more slowly than 
MgSiOa, and the gas density at which the condensation 
efficiency of Fe3C>4 grains becomes 0.5 is higher than for 
MgSiOs by a factor of ~5. This is because the volume of 
a monomer of Fe3C>4 grains (oc dp C3 o 4 o) is smaller than 
MgSiO-3 and thus the rate of grain growth is smaller (see 
Equation [2]) . This is also true for Mg2SiC>4 grains for 
which the volume per key element, Mg atom, is smaller 
than for MgSiOs. However, for these grain species, the 
critical metallicity is the same as for MgSiOs grains since 
these grain species have similar mean opacities (Figure 
[TJ) and thus similar cooling efficiencies through thermal 
emission. 

On the other hand, Si grains have mean opacity smaller 
than MgSiOs by orders of magnitude (Figure Q]). Thus, 
even if the condensation efficiency reaches unity, the gas 
cloud cannot cool down enough to meet the fragmen- 
tation conditions even for metallicities Z ~ 10~ 2 Zq. 
The growth of Fe grains results in larger critical abun- 
dance than MgSiOs because the species also have smaller 
opacity than MgSiOs. Hence, these grain species are un- 
likely to affect heavily the evolution of clouds. Indeed, 
Si and Fe grains could be rapidly oxidized into silicate 
and magnetite in a gas cloud with oxygen-rich compo- 
sition. Supernova models (N. Tominaga et al. 2013, 
in preparation) a nd observations in the Galactic halo 
(|Yong et al.l I2012D have suggested that the abundances 
of oxygen in metal-poor stars are usually Ao/Ax > 10 
where X = Mg, Si, and Fe. 

It is considered that C grains might have also large 
effects on the evolution of collapsing clouds. Yet, it is 
uncertain whether C grains can grow sufficiently in the 
collapsing gas because a large fraction of carbon atoms 
are expected to be locked in CO molecules and are no t 
available for the growth of C grains (jNozawa et al.ll2012l ) . 
Even if CO molecules are accreted onto dust nuclei, the 
CO ice mantle would be easily evaporated because of its 
low sublimation temperature. Nevertheless, if the initial 
condensation factor of C grains is high enough, the cool- 
ing of the gas by carbon dust thermal emission may have 
an impact on the cloud evolution. 

We have examined the effect of growth of an in- 
dividual dust species on collapsing gas clouds sep- 
arately. In a low-metallicity gas, multiple dust 
species would grow simultaneously, and thus their com- 
bined effects could reduce further the critical metallic- 
ity. Dust formation calculations i n primordial super- 
novae (e.g.. iTodini fc Ferrarall200U iNozawa et al.l 120031: 
ISchneider et al.ll2006HBianchi fe Schneiderll2007l f showed 
that various grains species can be produced. It is im- 
portant to study how the growth of multiple species of 
dust affects the thermal evolution of collapsing gas clouds 
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based on the realistic supernova dust models. 

Finally, we have set the sticking probability s» = 1 
in the calculations. As seen from Equation @, a lower 
sticking probability leads to a lower growth rate of dust. 
We perform another calculations, setting Si = 0.1. As 
a result, the fragmentation fails for Z — 1CP 5 Z & and 
r MgSiO 3 ,0 = 0.01 /zm. However, even for Sj = 0.1, grains 
with rM g sio 3 ,o = 0.001 /im rapidly grow to cause the 
cloud fragmentation. We have considered dust grains 
with a single initial size, but dust grains that were formed 
and processed in supernovae may have size distribution. 
Since smaller grains grow more rapidly, the lower limit of 
size distribution and the initial mass fraction of the small 



grains could be crucial quantities for cloud fragmenta- 
tion. We will investigate these effects comprehensively 
in our future work. 
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